Fermion- and Spin-Counting in Strongly Correlated Systems 
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We apply the atom counting theory to strongly correlated Fermi systems and spin models, which 
can be realized with ultracold atoms. The counting distributions are typically sub-Poissonian and 
remain smooth at quantum phase transitions, but their moments exhibit critical behavior, and 
characterize quantum statistical properties of the system. Moreover, more detailed characterizations 
are obtained with experimentally feasible spatially resolved counting distributions. 



I. INTRODUCTION 



A. Particle- and spin-counting 



Particle-wave duality is one of the most spectacular, 
and at the same time intriguing phenomena of quantum 
mechanics. Nevertheless, careful counting of particles, 
such as photons, in a given quantum mechanical state 
allows to fully reconstruct the wave nature and coher- 
ence properties of the state. The formulation of photon- 
counting theory in the frame of quantum electrodynamics 
[l| gave birth to modern quantum optics. Recent progress 
in physics of ultracold atoms made possible to develop 
and apply techniques of single atom counting to various 
systems. Since the pioneering experiments of Shimizu 0, 
spectacular measurements of Hanbury Brown - Twiss ef- 
fect for bosons [sj , and fermions 4] have been performed 
with ultra-cold meta-stable Helium atoms. Esslinger's 
group employed cavity quantum electrodynamics tech- 
niques to measure the pair correlation function in an 
atom laser beam outgoing from a trapped Bose conden- 
sate . These new detection methods allow in principle 
to measure full atom-counting distributions with spatial 
resolution (by counting only atoms in a certain spatial 
region), and provide novel efficient ways of detection of 
strongly correlated systems [6] . 



Equally spectacular progress has been achieved in 
spin-counting, or in other words, measurements of total 
atomic spin for atoms with spin, or pseudo-spin degree of 
freedom. The idea of quantum non-demolition polariza- 
tion spectroscopy (QNDPS), has been demonstrated in 
Ref. [3]. It employs the quantum Faraday effect: polar- 
ized light beam passed through the atomic sample, un- 
dergoes polarization rotation. Atomic fluctuations leave 
an imprint on the quantum fluctuations of the light, and 
vice versa. This idea was recently extended to ultra-cold 
spinor gases [sj, where it can be used to detect, manip- 
ulate, and even engineer various states of such systems. 
Amazingly, this method allows also for a spatial resolu- 
tion (when standing laser beams are employed) [^. 



B. Main results 

In this paper, we show how the atom counting tech- 
niques can be used to detect properties of strongly corre- 
lated systems. We concentrate, in particular, on the case 
of fermion and/or spin counting in one-dimensional (ID) 
optical lattices, that are equivalent, via Jordan- Wigner 
transformation 'l^, to ID spin chains. The problem of 
spin counting for a local block of spins in the ID Ising 
model in a transverse field has been considered in a beau- 
tiful work of Demler's group 11 ll . Our paper is in a 
sense complementary to Ref. |ll| . First, we consider 
not only on the Ising model, but on the whole family 
of asymmetric XY models, characterized by the asym- 
metry parameter 7, in the transverse field h. Second, 
employing ideas of Ref. we calculate not only the 
counting distribution for the total fermion number (total 
Z-spin component), but also for "effective" number, cor- 
responding to certain spatial Fourier components of the 
fermion density. While for the considered family of mod- 
els, counting distributions are always smooth, their cu- 
mulants exhibit critical behavior, evident even for small 
detection efficiencies. The distributions are always sub- 
Poissonian, but the sub-Poissonian character changes, as 
we sweep h from to 00. For small (large) 7, the h — 
(h = 00) distribution is always the narrower (broader) 
one. Through the paper, we use an elegant generaliza- 
tion of the photon-counting theory to fermions, derived 
by Cahill and Glauber within the formalism of Grass- 
mann variables [l^. For the cases we consider, we ob- 
tain analytic expressions for the counting distribution in 
terms of simple recursion relations. 

The paper is organized as follows. In Sec. |TT1 we briefly 
describe the models of the ID optical lattice that we con- 
sider, and the Jordan- Wigner transformation that can be 
used to diagonalize them. In the next section (Sec. IIIip . 
we derive the counting statistics of fermions in the sys- 
tems described by these models; in particular, we discuss 
them for the Ising model (Subsec. IIII Dp . and more gen- 
erally for the asymmetric XY model (Subsec. IIIIEp . In 
Subsec. IIII Fl we consider the means and variances of 
the counting distributions: We derive recurrence rela- 
tions that allow for easy calculation of these moments 
for an arbitrary number of particles in the system. We 
discuss also the generalization of our method to the case 
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of Fourier components of the total spin in Subsec. IIII HI 
Finally, we summarize our results in Sec. IIVI 



II. FERMI GAS IN AN ID OPTICAL LATTICE 



one should additionally introduce tunneling assisted with 
a laser or microwave induced double spin flip. For this 
aim, one should make use of the resonance between on- 
site two atom "up-up" and "down-down" states, without 
disturbing "up-down" configurations. 



A. ID Fermi gases 

Let us consider a family of models describing an one- 
dimensional Fermi gas in an optical lattice, described by 
the Hamiltonian 



j=0 



h.c. — 2gc'jCj 



(1) 

where J/2 is the energy associated to fermion tunneling, 
g = h/J, and N is the number of sites. One way to 
realize such Hamiltonian with ultracold atoms is to use 
a Fermi-Bose mixture in the strong coupling limit. In 
this limit, the low energy physics is well described by 
fermionic composites theory |13i | , in which fermions form 
composite objects with 0, 1, . . . bosons, or bosonic holes 
repectively. The fermionic composites undergo tunnel- 
ing and interact via nearest neighbor interactions, which 
may be repulsive or attractive, weak or strong, depend- 
ing on the original parameters of the system, such as 
scattering lengths, etc. In the case of weak attractive 
interactions, the system undergoes, at zero temperature, 
a transition into a "p-wave" superfluid, described well 
by the Bardeen-Cooper-Schrieffer (BCS) theory, corre- 
sponding exactly to the Hamiltonian ([1]). 



B. ID spin chains 

Using Jordan- Wigner transformation [l3|, one can 
transform the Hamiltonian ([1]) into the one of a ID asym- 
metric XY spin chain in the transverse magnetic field h, 



N-l 



+1 



(2) 

-'J I'^l ^""^ ^Pi"^ ^/'^ operators at site j, pro- 
portional to Pauli matrices. The special cases 7 = (i.e. 
the so called symmetric XY, or XX limit) and 7 — ±1 
can be realized with single species bosons in the hard 
core (i.e. strongly repulsive) bosons limit or in 

a chain of double well sites filled with bosons interacting 
via weak dipolar forces ^15,] , respectively. In general, one 
should use a two component Bose-Bose and Fermi-Fermi 
mixture, which, in the strong coupling limit, and in the 
Mott insulator state with one atom per site, is described 
by an asymmetric (XXZ) Heisenberg model (cf. [lH) in 
the Z oriented field. By appropriate tuning of the scat- 
tering lengths via Feshbach resonances, one can set the 
Sjj^^S^ couphng to zero, i.e. achieve the XX model in the 
transverse field. In order to introduce the asymmetry 7, 



C. Jordan- Wigner transformation 

As it is well known, Jordan- Wigner transformation 
works for open chains, and in particular for an infinite 
chain. We will nevertheless assume periodic boundary 
conditions to solve the fermion model ([T]) using Fourier 
and Bogoliubov transformations (see e.g. [13 )• For large 
N , such precedure gives the right leading behaviour. We 
define Fourier transformed operators as 



N-l 



4 = E exp(-«j«'/c)aj'. 



and 



fe=0 



N-l 



Cj = ^ exp(ij$fe)afe. 



(3) 



(4) 



fe=0 



where $fe = 2'Kk/N . We perform then the Bogoliubov 
transforms 

ak ^ Ukdk - ivkd}j^ a\ ^ Ukd\ ivkdN~k, (5) 
where u^, Vk are real numbers satisfying 



ul+vl = 1, 

UN-k = Uk and v^^k 



-Vk, 



so that we can write 



Uk = cos — , and Vk = sin — . 
2 2 



When 



tan 6* = 



7 sin $fc 
cos<I>fe - g' 



(6) 



(7) 



(8) 



the Hamiltonian reduces then to the noninteracting 
fermions Hamiltonian, 



= - XI ^k4dk, 



k=0 



with 



Ck = 2 J (cos $fc - 5)2 + 72 sin2 ^k 



(9) 



(10) 



The ground state is thus the vacuum of the dk operators. 
For 7 > the spectrum is everywhere gapped, except at 
the critical point gc — ^- For 7 = 0, dk's coincide with 
a^'s or a^'s, and the ground state is a Fermi sea. For 
— 1 < 5 < 1 the spectrum is then gapless and the system 
critical. Note that the number of original fermions Nf = 

^^Q^ clci, as well as the the total Z-component of the 



spin, = 
except at 7 



Nf 



1/2 are not conserved. 
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III. FERMION COUNTING STATISTICS 
A. Fermion counting distributions 

Let us now turn to counting procedures. For the case of 
fermions, one should think about the analogue approach 
as one used in the experiments on metastable Helium. 
For spins, one could use directly QNDPS to measure the 
distribution of Sz, or even its spatially resolved version 
[9] . An alternative way would be to switch off the Hamil- 
tonian ^ (by switching off lasers), and induce sponta- 
neous Raman transition from the state "up" to some side 
level. Counting of spontaneously emitted photons would 
correspond then to counting of "up" spins" 

Mathematically, as known for photons and gener- 
alized by Cahill and Glauber for fermions , the prob- 
ability of detecting m photons in a given interval of time 
can be expressed as the mth derivative with respect to a 
parameter A of the generating function Q(A) as 



p{m) 



(-ir 



-Q 



(11) 



where Q(A) is the expectation value of a normally or- 
dered exponential Q(A) = Tr(p : e~^-^ :). The oper- 
ator I is a space-time integral of the product of the 
positive-frequency and negative-frequency parts of the 
quantum fields describing particles to be counted. The 
mean values of normally ordered products can be calcu- 
lated in a particularly convenient and elegant way using 
the Grassmann variables formalism, introduced in p^ . 
In the case of counting the total number of particles, we 

have I = kJ^j^o ^Kj2j=o ^'j'^i = '^Efc^o «fc«fc > 
where k,— \ — exp(— ^t), while < C < 1 is the detector 
efficiency, and t is the exposure time. For the spatially re- 
solved QNDPS, X = i^^Y^=Q cos(kirj), where k^, is 
the wave vector of the standing wave used for detection, 
and Vj is the position of the j-th site. 

For counting the total number of particles, the gener- 
ating fmiction Q(A) can be written as 



Q(A) = Tr(p : e^^^^ ^"=0' :). 



(12) 



The operators a^a/j commute for different fc, so that the 
expression for Q can be rewritten as 



AT-l 



Q(A) = Tr(p: [] (e^^^'^^'^'' ) :) 

Af-l 

= Tr(/3 : (1 — AKa|afc + A^K^a|afea[.afc -|- .. 

fc=0 
N-l 

= Tr(p ]^ (1 - \Ka\au)) 



)■■) 



k=0 
N/2 



Tr(p JJ(1 - AKa^afe)(l - XKalj_^aN-k)), 



as : al.akal.ak := al.al.akak = 0, etc. 

The terms a^at- and a^^_i^aN-k can then be expressed 
in terms of the d fermions: 

alak = {ukd\ + ivkdN^k){ukdk - ivkdl^^k)' 
ajv-feOjv-fc = ("fcc?]v_fc - ivkdk){ukdN-k + ivkd\). 

B. Generating function for the ground state 

We consider the counting statistics of the c fermions in 
the ground state of the Hamiltonian, i.e. in the vacuum 
state of d fermions. 

The trace in the generating function can be now easily 
calculated by the formalism of Grassmann variables . 
The P representation for the density operator p is 



(13) 



p= d dp{d)\d){d\, 



where |a) are the fermionic coherent states, as defined 
in ^13]. Using the P representation, the mean values 
of normally ordered products of d-fermions can then be 
calculated as 

Tr(p4"dr)= [ d^dP{d){d\di-dT\d) 

d^dP{d)al''aT, (14) 



where the are Grassmann variables, and are defined 
by the eigen-equation d\ai) = ai\ai). For the vacuum 
state of the d-fermions, 



p= |0...0)(0...0|, 

the P-function is given by 



(15) 



Pia) 



J d^Cexp {Y.i^^^: - = S{d). (16) 

Evaluating Eq. (fT4|) using Eq. (fT6|) . we get the relations 

TTipdi-dT) = f d'dHiairarSid) = 0, (17) 
•' i 

and 

Tr(p) = / d^dS{d) = 1. (18) 



k=l 



The relevant remaining terms in the product (1 — 

XKa\ak){i~ XKa''j^_i^aN-k) in the generating function are 
thus 

9 t 9 t 

1 - XnVkdN-kd^M^k ~ Xnvkdkdl 
+X^K^vldN-kd^^_kdkd\ ^ X^K'^vlu'idN-kdkd\_kd\ 



4 



Elementary calculations using the relations p?|) and (fT8|) 
yield 



N/2 



(19) 



fc=i 



At this point it is convenient to introduce the distribution 
function p{m, M) of counting m particles for M pairs of 
modes. It is given by the same expression as before, but 
with the product in Eq. limited to M/2 terms. 



C. Counting statistics 

The counting distribution is calculated from the gen- 
erating function by the relation in Eq. (fTTj) . We use the 
generalized Leibniz rule, 



fc=i 



E 



N 



where the generalized Newton's symbol is given by 



rii, n2, njv/ ni!n2!...nAr! ' 

to derive a recurrence relation, to calculate the distri- 
bution for (M -I- 1) modes, given the distribution for M 
modes. 

The distribution function p{ni, M) for M modes is 
given by 



p(m, M) 
(-1)™ 



{-ly 



E 



-Q 

■ Ia=i 

M/2 . 



h\l2\...lN\ 



(20) 



where the summations run over li^ ... ,1m such that li 



■ ■ ■ + Im = 'TT-j where = 0, 1, or 2, for j = 1, . 
We can now derive the recursive relation 



2 

E^ 

1=0 



p{m,M+l) = V^'Pip(m - i,M) 



(21) 



where 



Vo = 1 - 2kw^,^^i -I- K^w^.^+i, 



7^2 = l-T'o-T'i 



(22) 



are the probabilities of detecting 0,1, or 2 particles in 
the modes M + 1 and N ~ M ~ 1. Therefore, start- 
ing from p(0, 1) = 1 — 2kv1 + K^vf, p{l, 1) = Anvf and 
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FIG. 1: Counting statistics of the transverse Ising model. 
The horizontal axis is (m — rh)/N + 1. The vertical axis 
is of the corresponding probability. The curve with purple 
squares is for h/ J = 0.01, while the one with green circles is 
for h/J = 10. The QPT of this model is at h/J = 1. Both 
the distributions are sub-Poissonian. However, the counting 
distribution becomes much narrower in the case when h/ J > 1 
than the situation when h/J < 1. In this case, we have taken 
the efficiency k as 0.9. 



we can use the recurrence relation (|2T]) to 



pi2,l)^Kvl 

calculate the counting distribution for an arbitrary num- 
ber of modes. 

Let us turn now to our results and discuss the counting 
statistics for different values of 7. In the figures that we 
plot below (except in Fig. [6] in Subsec. IIII Gp . we choose 
a value of the total number of modes, N, such that the 
corresponding quantities (distribution, mean, variance, 
etc.) have already converged. In the cases that we con- 
sider, such convergence occurs for N « 300. 



D. Transverse Ising model 



The counting distributions for the transverse Ising 
model (transverse XY model with 7 = 1) for two ex- 
emplary values of the field parameter g = h/ J are shown 
in Fig. [TJ The Ising model has a quantum phase transti- 
tion at g = 1 [lol | , and one exemplary value of g is chosen 
below the QPT, and the other above it. The difference 
in behavior is clearly seen, (to and var denote the mean 
and variance of the distribution, respectively.) Below, it 
will be more clearly revealed by looking at the mean and 
the variance of the distribution. 
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which can be calculated from the following recurrences: 



! 1 



/I 
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FIG. 2: Fermion counting distribution as a function of (m — 
fn)/N + 1 (horizontal axis) for k — 0.9, and for the indicated 
values of 7. Purple squares correspond to /i/J — > 0, while 
green circles to h/J — > 00. The transition anisotropy is here 
at 7 f« 0.1. 



E. Transverse XY model: "Transition anisotropy" 

In Fig. [21 we plot counting distributions as a function 
of (m — rn)/N + 1 for four values of 7, for a fixed value 
of the efhciency k — 0.9, and for two extreme values 
of 3: g — > and g 00. Note that all the distribu- 
tions presented in Fig. [5] are smooth and their widths 
(~ ^/var/N) are of order of 0.01. Since, as we argue be- 
low, m ~ kN, all the distributions are sub- Poissonian, 
because var < in, despite the finite detection efficiency. 
For 7 — !■ 0, the distribution for g — )■ is narrower than 
that for g 00. This tendency is inverted in the Ising 
model, when the distribution for .g — > has a larger vari- 
ance than the one for g ^ 00. At, what we call, transi- 
tion anisotropy 7 ~ 0.1, the distributions for 5 — s- and 
g ^ 00 practically coincide. 

This transition anisotropy depends on the efficiency /t, 
and it moves to 7 — > 0, as k — > 1. This indicates that 
the probability distribution of counting can distinguish 
the two universality classes (the XX, with 7 = 0, and 
the Ising, with 7 > 0) among the XY models on a chain. 
In the limit of k — > 1—, only the model with 7-^0 has 
lower variance for (7 — > as compared to g ^ 00, while 
all the other XY models (with 7 7^ 0) have the opposite 
behavior. 



F. Recurrence relations for mean and variance 

In order to understand the properties of counting dis- 
tributions better, we look at the mean and variance. 



2 

inn ^ 



varM+i = "T-M+i - "T-M+i 



varM + 4k t;M+i(l 



1) 



(23) 



(24) 



Since mT and vari can be trivially calculated, the mean 
and variance can be obtained by these relations for an 
arbitrary number of modes. The recurrences imply that 
the mean vTm < kTV; we find typical value of ttTtv in- 
deed of order of kN . On the other hand, the variance 
war^v < K^N . Both quantities show singular behavior in 
the thermodynamical limit at criticality. In particular, 
for the transverse Ising model (7 = 1), near the critical 
point g — gc = 1, the mean m can be written in terms 
of elliptic integrals of first and second kind, and can be 
expressed as [l^ 

m w ~-^{g - gc) In Iff - 9c\ - 

it: it 



so that 



dfn/dg w —(In I5 — 5c| + l)/27r. 



Since all models with 7 7^ belong to the same univer- 
sality class, they all present the same singular behavior 
This is contrasted with the case of XX model, which 
belongs to a different universality class. The singular be- 
havior is clearly seen in the plots of m/N and var/N 
obtained for finite N ~ 300 and ideal k = 1 (see Fig. 
[21). For finite values of 7, the variance shows a jump in 
the first derivative, while the first derivative of the mean 
tends to "infinity" at gc- This behavior is better seen, 
when one plots directly the derivatives of m and var (see 
Fig. U). This behavior changes drastically as 7 0. 
The variance tends then to zero (in the symmetric XX 
model the particle number is conserved), and the mean 
has a diverging derivative for g < gc, and is constant for 
g > gc- Amazingly, although finite detector efhciency 
obviously smoothes out the curves, the signatures of the 
singularities are clearly visible even for k = 0.5 (see Fig. 

A clear change of behavior of the curves is visible 
even at k — O.l! Note, that in all considered cases so 
far, the variance var/N < rn/N, i.e. all distributions 
are sub-Poissonian. Note, however, that going from anti- 
ferromagnetic to the ferromagnetic case, does not affect 
the variance, but replaces m/N {1/2 — m/N)- In that 
case we do observe a transition from sub-Poissonian be- 
havior at small g < gt, to (weakly) super-Poissonian for 
g > gu the value of gt tend to gc from below as 7 ^ 0. 



G. Even versus odd splitting 

The Bogoliubov transformation used to solve the con- 
sidered models can be regarded as a "squeezing" or "pair- 
ing" transformation. The ground state that we investi- 
gated is analogous to BCS states of semiconductors, i.e. 



6 




0.8 
0.6 
0.4 
0.2 

0,0 -k« ««««»««»»»»€. 



r=o.oi 


0.4- 




7=0.5 








K=^ 




0.2- 









2 



1 

h/J 



0.4 



0.2 




2 




0.08 



0.04 



0.00 




FIG. 3: Mean m/N (blue squares) and variance var/N (red FIG. 5: Mean m/N (blue squares) and variance var/N (red 
circles) of the fermion counting distribution as a function of circles) of the fermion counting distribution as a function of 



g — h/J for At = 1, and indicated values of 7. 



g = h/ J for 7 = 1, and indicated values of k. 



ji 1.2- 7=0.5 ■ 

1.0- i':^ - k=l I 1.0- 




1 20 1 , ,, 2 1 2 



FIG. 4: The derivatives of the means and variances are plotted 
against the transverse field h/J (horizontal axis), for 7 = 0.01, 
7 = 0.5, and 7 = 1. Blue squares denote the derivatives 
of the means, while red circles denote the derivatives of the 
variances, in the respective cases. Also, k = 1. The QPTs of 
all the models at y = 1 are clearly visible. 




3000 3050 3100 31 50 3200 3250 3300 4750 4800 4850 4900 



FIG. 6: Even versus odd splitting for n = 0.999 in the Ising 
model. For N=1000 the probability distribution splits up, 
whereas for N=4000 there is no splitting. 



they involve fermion (Cooper-like) pairs. Thus, in the 
ideal case of k = 1, the counting distributions are ex- 
actly zero for odd numbers of particles. In practice, for 
finite values of A'' and k < 1, the distributions oscillate 
between larger values for even, and small for odd num- 
ber of counts. This behavior is very strongly affected by 
K < 1, since at finite efficiency, one may easily miss sin- 
gle atoms from the Cooper pairs, and obtain odd counts. 



In effect, for a given value of TV, the even-odd asymme- 
try is visible only for k close enough to 1. Similarly, the 
even-odd asymmetry is strongly affected by the finite size 
effects - for a given value of k < 1 it is visible only for N 
small enough (see Fig. [S]). 
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H. Counting spatial Fourier components of the 
fermion density 

Finally, let us point out that the methods proposed 
in allow for measurements of various kinds of Fourier 
components of the total spin; in terms of particle count- 
ing, these methods allow for instance to count particles 
in every second, every third site, etc. Our theory is easily 
generalized to such situations. 

In the case when we count every second c fermion, we 
have to express &2j^2j = '^tj'^'^j terms of the d fcrmions. 
As before, as a first step we do the Fourier transform: 

N-l 



N-l 

C2j = ^ exp(+2?j$fc)afe. 



k=0 



N/2-1 



The expression J2j=o 
be written as 



t _ 1 Y^«- 
'^2j'^2j — 2 2.^7=0 



N-l t 



■2j 



(25) 



C2j can thus 



N/2-1 
3=0 



1 1 - exp (47r«(/c - fc')) 

2 ^ r 



k,k' 



exp (47ri(A: - k')/N) 



alak', 



which is non vanishing for k 
Finally 



fc' = or Ifc - fc'l = ^- 



(26) 

N 
2 • 



(27) 



We can now calculate Q(A) as follows: 

JV/2-1 



Q(A) = Tr(/9 : Y\ e^^^''''"'="'"°"/2+'^^'''"'"'"''"^^+''-' :) 

k=0 
N/2-1 ^ 

= n 2^'^(^k+(^N/2+k)i(^k+aN/2+k)) 
N/A ^ 

= n ((^ ~ 2^'**-"'= "jV/2+fe)(afc + aN/2+k)) 
k=l 

X (1 - ^^'«(ajv-fc + ajv/2-A;)(aJv-fe + aw/2-fc))) ■ (28) 

After performing the Bogoliubov transform, and keep- 
ing the relevant terms for the vacuum state of the d 
fermions, the generating function Q is given by 



JV/4 



(29) 
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FIG. 7: Mean m/N (blue squares) and variance var/N (red 
circles) of the counting distribution of every second fermion 
as a function of g — h/J for 7 = 1, and indicated values of k. 



which is in the same form as in Eq. (|19p . with the prod- 
uct restricted, however, to — 1 terms. We then easily 
derive analogous recurrences as in the cases considered so 
far. Fig. [7] show the behavior of the mean and the vari- 
ance, when counting every second spin, in the transverse 
Ising model. Note that the traces of singular behavior at 
g — gc persist. What is perhaps more interesting is that 
the general behavior is more rich. In particular, there 
is a crossing from sub- to super-possoinian behavior at 
0.5. For 7 — !■ the point of crossing moves to zero, 
and the variance disappears. 



IV. SUMMARY 



Summarizing, we have formulated and applied fermion 
and spin counting theory to a family of one-dimensional 
strongly correlated systems that can be realized and de- 
tected with ultracold atoms. The counting distributions 
exhibit traces of singularities at criticality, that persist 
even at low detection efficiencies. They show various 
kinds of rich behavior, such as transitions from sub- to 
super-Poissonian character and even-odd oscillations. 
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